Jet modification in three dimensional fluid dynamics at next-to-leading twist 
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The modification of the single inclusive spectrum of high transverse momentum (pr) pions ema- 
nating from an ultra-relativistic heavy-ion collision is investigated. The deconfined sector is modelled 
using a full three dimensional (3-D) ideal fluid dynamics simulation. Energy loss of high px par- 
tons and the ensuing modiflcation of their fragmentation is calculated within perturbative QCD at 
next-to- leading twist, where the magnitude of the higher twist contribution is modulated by the 
entropy density extracted from the 3-D fluid dynamics simulation. The nuclear modification factor 
(Raa) for pions with a pr > 8 GeV as a function of centrality as well as with respect to the reaction 
plane is calculated. The magnitude of contributions to the differential Raa within small angular 
ranges, from various depths in the dense matter is extracted from the calculation and demonstrate 
the correlation of the length integrated density and the Raa from a given depth. The significance 
of the mixed and hadronic phase to the overall magnitude of energy loss are explored. 

PACS numbers: 12.38.Mh, ll.lO.Wx, 25.75.Dw 
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Current experiments at the Relativistic Heavy Ion Col- 
lider (RHIC) have established two fundamental proper- 
ties of the dense partonic matter produced in the col- 
lision of heavy ions The produced matter achieves 
local thermal equilibrium rapidly and is endowed with a 
very low viscosity; as a result, the dense partonic matter 
demonstrates almost ideal hydro-dynamic expansion 0. 
Secondly, the propagation of high transverse momentum 
partons, produced in the early hard scatterings, is sys- 
tematically impeded. The spectrum of high transverse 
momentum hadrons with Pt > 8 GcV shows a reduc- 
tion of almost a factor of 5 in central Au-Au collisions 
compared to that from the spectrum of such particles in 
elementary nucleon nucleon encounters enhanced by the 
number of expected binary collisions 

In most estimates from jet-quenching calculations 0, 
[1], there exists no real dynamical model of the 
medium Q; in the best of circumstances, a time depen- 
dent, three dimensional, functional form p(x,y, z,t) ^ 
Po{x, y, z)(to/t), with the initial density estimated from a 
woundcd-nuclcon model, is used. There have been exten- 
tions involving a 1-D hydrodynamical model , however 
such calculations cannot address the centrality depen- 
dence of the modification. The first attempt involving 
a realistic model of the medium 0] estimated the effects 
of 3-D expansion on the uncharged pion spectrum and 
the Raa- However, this approach treated the energy loss 
of jets in a rather simplified manner. Detailed 3-D sim- 
ulations of the medium have recently become available 
in Ref . 0] , where the various parameters of the simula- 
tion were fitted by comparison with experimental data 
in the soft sector. A study of the modification of jets 
in this medium, using the quenching weights of Ref. Q, 
was carried out in Ref. [loj . However, such calculations, 
along with others have tended to reformulate the effect of 
a 3-D expanding medium in terms of an effective length. 
In the current effort, the medium modified fragmenta- 
tion functions of high pT pions will be computed in the 



higher- twist energy loss formalism llj . The 3-D expand- 
ing plasma of Ref. 0] will furnish the space-time depen- 
dence of the two-gluon matrix element (related to the 
well known transport coefhcicnt g), which parameterizes 
the sensitivity of the jet-modification equations to the 
properties of the medium. The calculation of the contri- 
butions which lead to the Landau Pomeranchuck Migdal 
(LPM) interference will be carried out directly in such 
a medium without recourse to an effective length. In 
the following, the salient features of the 3-D fluid dy- 
namical simulations and the ensuing calculation of the 
nuclear modification factor at next-to-leading twist will 
be briefly reviewed. 

Ideal Relativistic Fluid Dynamics (RFD, see e.g. (l^ ) 
is the study of matter where the precise degrees of free- 
dom are unimportant to the evolution of the extensive 
variables of the system. It is ideally suited for the high- 
density strongly coupled phase of heavy-ion reactions at 
RHIC, but breaks down in the later, dilute, stages of the 
reaction when the mean free paths of the hadrons be- 
come large and flavor degrees of freedom are important. 
Relativistic Fluid Dynamics is based on the solution of 
the four equations of energy momentum conservation, 

d^T^'^ix) = 0, (1) 

where T^'^ is the energy momentum tensor, given by 
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Here, e, p, U and g^'' are the local energy density, pres- 
sure, four velocity and metric tensor, respectively at x. 
The relativistic hydrodynamic equation [Eq. ^] is solved 
numerically using baryon number ub conservation 



^^{nB{T,^l)U^') = 0, 



(3) 



as a constraint and closing the resulting set of partial 
differential equations by specifying an equation of state 
(EoS): e ^{p)- In fhs ideal fluid approximation (i.e. 
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neglecting off-equilibrium effects), once the initial condi- 
tions for the calculation have been fixed, the EoS is the 
only input to the equations of motion and relates directly 
to the properties of the matter under consideration. As 
a result, even striking phenomena occurring within sys- 
tems in local equilibrium, such as phase transitions may 
be easily incorporated into the bulk dynamics of the full 
system. Ideally, either the initial conditions or the EoS 
should be determined beforehand by an ab-initio calcu- 
lation (e.g. the EoS may be obtained from a lattice-QCD 
calculation), in which fit to the data would allow 

for the determination of the other quantity. 

We assume that hydrodynamic expansion starts at 
To = 0.6 fm. Initial energy density and baryon number 
density are parameterized by 



nB{x,y,r]) 



M{x,y;b)H{r,), 
,axW^(x,y;&)i?(7?), 



(4) 



where b and Cmax (JT-Smax) are the impact parame- 
ter and the maximum value of energy density (baryon 
number density), respectively. W{x,y;b) is given by 
a combination of wounded nuclear model and binary 
collision model [l^ and H{ri) is given by H{r]) = 
exp [— (|?7| — ?7o)^/2cr^ • 9{\ri\ — 779)] ■ RFD has been very 
successful in describing soft matter properties at RHIC, 
especially collective flow effects and particle spectra 
Il4l |. All parameters of our hydrodynamic evolution 9[ 
have been fixed by a fit to the soft sector (elliptic flow, 
pseudo-rapidity distributions and low-px single particle 
spectra), therefore providing us with a fully determined 
medium evolution for the hard probes to propagate 
through. While different descriptions for the hadronic 
phase were presented in Ref. 0], in this calculation, the 
entire evolution of both the deconfined and hadronic 
phase will be assumed to be describable by the RFD sim- 
ulation. 

Assuming a factorization of initial and final state ef- 
fects, the differential cross-section for the production of a 
high pt hadron at midrapidity from the impact of two nu- 
clei A and B is obtained as a convolution of initial nuclear 
structure functions, G^(xa, Q^), G^jxb, Q^) (shadowing 
functions are taken from Ref. [ia|)j a hard cross section to 
produce a high px parton (da / dt) and a medium mod- 
ified fragmentation function. The entire effect of final 
state jet modification is encoded in the medium modified 
fragmentation function D(z,Q^) expressed as the sum 
of the leading twist vacuum fragmentation function [l6| 
and a correction from re-scattering of the struck quark 
in the medium i.e., D — D + AD In the coUinear 

limit, the modification is computed by isolating correc- 
tions, suppressed by powers of Q^, which are enhanced 
by the length of the medium 17[. At next-to- leading 
twist, the correction for the fragmentation of a quark has 
the expression (generalized from deep-inelastic scattering 

(Dis) P[i3), 
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X llN,Gtixa)Gi{x,)^ 
dt 



+ V.C. (5) 



In the above equation, l± is the transverse momentum 
of the radiated gluon (quark) which leaves a momen- 
tum fraction y in the quark (gluon) denoted as parton i 
[Pq^i{y) is the splitting function for this process] which 
then fragments leading to the detected hadron. The fac- 
tors tA(f),tB(f) are the thickness functions of nuclei A 
and B at the transverse location f = {x,y). The v.c. 
refers to virtual corrections. The factor T*'^, which has 
its origin in higher twist matrix elements of the quark, 
may be expressed as. 



T^" = d'bd'rtA{r + b/2)tB{f- b/2)G^{xa)G^ {x^) 



da 
dt 



X / dCxgPg{xg,nC + r){2-2coa{r]LO)- (6) 
Jo 

In the above equation, the factor rjL = l']_/[2pTy{l — y)]. 
Where, px represents the transverse momentum of the 
produced parent jet. 

The jet direction (n) is chosen as the x-axis in the 
transverse (a;, y) plane (z-direction is set by the beam 
line); the angle of the reaction plane vector b is measured 
with respect to this direction. The distance travelled by 
the jet in this direction prior to scattering off a gluon 
is denoted as (. The gluon's forward (towards the jet) 
momentum fraction in denoted as Xg. In what follows, 
the nuclear modification of the single inclusive spectrum, 
both integrated and as a function of the angle with re- 
spect to the reaction plane will be computed. The angle 
integrated modification, the Raa, is defined as. 



Raa 



dyd'-'pT 



, \ daPP{pT,y) 
'mim'Jmaxj dyd^pr 



TAA{b 



(7) 



where, TAA{i>mim b,nax) is related to the mean number of 
binary nucleon nucleon encounters in the range of impact 
parameter chosen ((iVf,m)), the geometric A- A cross sec- 
tion [aQfg{bmin, bmax)] and the inelastic p-p cross section 

(cr™) as TAA{bmin,bmax) ^ (Nbin) a^^^{brmn , bmax) / a^^p . 

The phenomenological input that is required to un- 
derstand the variation of the nuclear modification factor 
with centrality or pT is the space-time dependence of 
the gluon density p{x,y,z,T) = XgPg(xg,x,y, z,t). This 
density is directly proportional to the transport coeffi- 
cient qji (for a jet parton in representation R) p^. [20|. 
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In the remainder, the z-coordinate will be set to zero as 
the interest will remain on midrapidity observables. The 
gluon density at a given location depends on the phase 
of matter at that location and is in general assumed to 
depend on the number of degrees of freedom prevalent in 
the excited matter at that location. The exact degrees of 
freedom in the dcconfincd phase are as yet unknown [2l| . 
The partonic density as seen by the jet within each de- 
gree of freedom is also unclear. For this first attempt at 
higher twist energy loss in a hydrodynamically expanding 
medium, we invoke the very simple ansatz. 



q{x,y,z,T) ^ qo -3 (8) 

^0 

X [R{x, y, z, t) + chg {1 - R{x, y, z, t)}] , 

where, T(x,y,z,T), ^i_(x,y,z,T) and R(x,y,z,T) repre- 
sent the temperature, flow transverse to the jet and the 
volume fraction in the plasma phase at the space-time 
point x,y,z,T. It is this information that is extracted 
from the RFD simulation. The factors QojTq represent 
the maximum q and temperature achieved in the simula- 
tion; in this particular version of RFD, Tq = 0.405 GeV 
and qo is a fit parameter adjusted to fit one data point of 
the Raa, at one centrality. The factor chg accounts for 
the fact that the partonic density per degree of freedom 
in the hadron gas phase may be different from that in 
the QGP phase and plots with different choices of chg 
will be presented. It should be pointed out that contrary 
to the (pre-)hadronic scattering picture of Ref. [l^l, the 
energy loss of jets in the hadronic phase is considered 
here as a completely partonic process where the jets 
fragment into hadrons outside the medium. 
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FIG. 1: (Color online) Raa (with qo ~ 1.3 GeVVfm and 
Chg = 1) in Au-Au collisions at 0-5% (blue line) and 20-30% 
(green line) centrality compared to data from PHENIX [2^ . 

The above prescription of assuming q to be a linear 
function of may be referred to as the thermal prescrip- 
tion. In the energy density prescription used in Ref. [l^l , 
one assumes g to be a linear function of e^/^. Alterna- 
tively, one may also assume that g is a linear function 
of the entropy density s. Relationships may be drawn 



between these various schemes with a knowledge of the 
EoS. For instance, in a plasma phase which is populated 
by massless partons, the entropy density is a linear func- 
tion of and the two prescriptions are identical. This 
breaks down in the hadronic phase as the entropy density 
is no longer a simple linear function of . Due to the va- 
riety of unknowns regarding the q in an excited hadronic 
resonance and as a result, its parametric dependence in 
the hadron resonance phase, we persist with the simple 
prescription of Eq. ([5]). 

As mentioned above, the overall fit parameter q^ is 
tuned to fit one experimental data point, at one central- 
ity and pt- For the current effort, the fit parameter is set 
by requiring that the Raa at pr = 10 GeV in the most 
central event (0-5% centrality) is 0.2. For a chg ~ 1 this 
corresponds to a ^0 — 1-3 GeV^/fm for gluon jets. With 
the value of go and chg fixed, the variation of Raa as 
a function of pt and centrality of the collision are pre- 
dictions. These are presented in Fig. [T] for two central- 
itics (0-5% and 20-30%) in comparison with data from 
PHENIX [2^. An important caveat to the estimates of 
q are the estimates of the plasma temperature from the 
very early phase of the evolution. Hard jets are produced 
very early (t < 0.1 fm/c) in the evolution of the system. 
However, fits to the spectrum of soft particles indicate 
that a locally equilibrated matter does not arise before 
0.6 fm/c. This is the starting point of the RED sim- 
ulations used. In the regime between 0.1 to 0.6 fm/c, 
estimates of q arc open to speculation. In this first ef- 
fort, the initial q is assumed to be held constant from the 
time the jets are produced to 0.6 fm/c. In future efforts, 
variations of this form will be attempted. 



0.5 



0.4 — 



ee2 



i 0.3 



0.2 




0.1 — 



■ ■ 




'=HG= 


= 1,(0-5%) 




'=HG= 


= 1(20-30%) ■ 




• 


'=HG= 


= 10,(0-5%) , 






'=HG= 


=10,(20-30%) 


♦ ■ 


♦ 


'=HG= 


=0,(0-5%) 






'^HG^ 


=0,(20-30%) ' 



10 



12 



14 

Pt 



16 



20 



FIG. 2: (Color online) Raa as a function of centrality for 
different choices of chg ~ inalT'^ in the hadronic phase. 
See text for details. 

In the RFD simulations used, an equation of state with 
a first order phase transition is included. As a result, 
such simulations include a prominent hadronic and mixed 
phase. In Fig.[TJ, the factor chg = 1 meaning that the q 
per unit of in the deconfined phase is similar to that in 
the hadron gas phase. From the equation of state used, 
one may discern the approximate and simple relation: 
sqgp ~ ^shgi between the entropy densities in the de- 
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confined phase {sqgp) and the hadronic phase {shg)i in 
the vicmity of the phase transition. This impHes that 
a chg = Ij places the ratio q/s, i.e., the quenching 
strength per unit of entropy, in the hadronic phase to 
be five times larger than that in the QGP phase. Due 
to the error bars in the data, a range of values of chg is 
allowed; however, certain extremes may be ruled out. 
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FIG. 3: (Color online) Raa as a function of the reaction plane 
for Au-Au collisions at 20 — 30% centrality. 

As an illustration, the Raa for the same centralities 
as in Fig. [1] are plotted for two other cases of c^g = 
and Chg = 10 in Fig. [21 In each case, the overall normal- 
ization factor go is re-adjusted to the standard condition 
that the Raa for a pion pT = 10 GeV and at a cen- 
trality of 0-5% be 0.2. This is the reason that all three 
plots for the three different choices of chg at 0-5% cen- 
trality essentially lie on top of one another. However, 
the Raa for the semi-central case (20-30%) shows no- 
ticeable variations. The red squares correspond to the 
case of Fig. [TJ The case Chg = 0, represented by the 
green diamonds, represents the case of no quenching in 
the hadronic phase. Due to the proportionately larger 
presence of the hadronic phase in more peripheral colli- 
sions, this case, understandably demonstrates less mod- 
ification [2^. The other extreme is with a chg ~ 10, 
where the q in units of is an order of magnitude larger 
in the hadronic phase compared to the dcconfined phase. 
As would be expected, the results show an opposite trend: 
the centrality dependence of the Raa is under-predicted. 

In semi-central events, the average path length trav- 
elled by jets in the reaction plane are different from the 
path lengths travelled out of plane. This should lead 
to the occurrence of an azimuthally dependent nuclear 
modification factor. Full three dimensional RFD sim- 
ulations capable of reproducing the observed azimuthal 
asymmetry in the low spectrum require that the ini- 
tial gluon density profile display a considerable azimuthal 
eccentricity in transverse space. The hard scatterings, 
which lead to the production of back-to-back jets tend 
to produce them isotropically. As a result, a sizeable 
azimuthal asymmetry in the nuclear modification fac- 
tor demonstrates the sensitivity of this observable to the 
early density profile of the QGP. As an illustration of 



this effect, the nuclear modification factor at an angle 
0° < < 90° with respect to the reaction plane, in an- 
gular bins of 15° JJ^j, is plotted in Fig.O One notes, that 
the Raa in plane shows the least modification, which 
grows as the angular bin is turned away from the reac- 
tion plane. As a result, it may be argued that the initial 
profile of the produced matter plays an important role in 
the determination of the final nuclear modification fac- 
tor. The rate of change of the Raa as a function of 
shows the interesting pattern that is slow at the terminal 
points and maximal in the middle of the angular range 
(i.e., ~ 45°). The results obtained from this analysis 
(and especially the symmetric behavior of the modifica- 
tion as a function of angle with respect to the reaction 
plane) show an identical pattern as that of the exper- 
imental measurements of the Raa versus the reaction 
plane at lower pT, as reported in Ref. [23|. The shallow 
rise in the Raa as a function of the pt of the detected 
hadron is demonstrated by the nuclear modification fac- 
tor in each angular bin: this property is thus independent 
of the density of the matter probed by the jet. 

NjjN^, are arbitratry normalizations 
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FIG. 4: (Color online) Raa as a function of the location along 
the impact parameter from longitudinal strips perpendicular 
to h. Also plotted in the number of initial binary nucleon- 
nucleon collisions and the length integral of the time depen- 
dent density weighted with the number of binary collisions 
(J dlp{x, y, z, t)) arbitrarily normalized to fit on the same plot. 

The existence of a calculation scheme for the modifi- 
cation of jets in a three dimensional dynamically evolv- 
ing medium allows for a more differential analysis of jet 
quenching in medium. An an example, we compute the 
Raa as a function of location in the reaction plane of the 
20-30% centrality collisions. The focus is on the Raa 
in the angular bin of 0-15° with respect to the reaction 
plane, this essentially corresponds to a 30° opening angle 
around the impact parameter b. The tansverse plane in 
then divided into ten longitudinal strips of width 1 fm 
perpendicular to b (as illustrated in the inset in Fig. d]). 
The Raa for jets originating in a given strip is then plot- 
ted (red dot-dashed line) as a function of the location 
of the mid-point of the strip in the reaction plane. Also 
plotted is the initial number of binary nuclcon-nucleon 
scatterings in each of the strips (black solid line). Plotted 
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along with this is the averaged length integrated density 
p [which is directly proportional to g, see Ref. [l^ H^] ] 
as experienced by a jet originating in each strip (green 
dashed line). This is defined as the following quantity 



{pL) 



(Pb(frtAir + b/2)tBif- 6/2) 



X / d(f>f,dCpir + Cn,T = C) 



(9) 



d^bd^rtA{r + b/2)tB{r- b/2) 



In the above equation, the density p is evaluated at the 
location r + and at time where r is the location of 
the jet vertex, C is the time elapsed and n is the direction 
of propagation of the jet. The range of (j>n is from —15° to 
15°, while r ranges within the specified strip. Thus, (pL) 
represents the length integrated density as experienced 
by a jet on average in the given strip. As may be noted 
from Fig. HI the Raa hi a given strip is anticorrelatcd 
with the mean integrated density experienced by the jet 
on its way out of the medium. The Raa tends to show 
a small rise in the strip farthest from exit (note that this 
is correlated with a drop in (pL)), this is due to the fact 
that a considerable fraction of jets that originate on this 
strip assume tangential paths of exit that do not pass 
through the center of the dense system produced. As 
one moves towards the middle of the impact parameter. 



a large fraction of jets, due to the angular cut, are forced 
to pass through the densest part of the matter and as a 
result, the Raa first drops and then rises. 

In this letter, higher twist calculations of energy loss 
have been extended to include realistic medium evolu- 
tion as stipulated by a 3-D hydrodynamical simulation. 
The nuclear modification factor at different centralitics 
as a function of the pr of the detected hadron has been 
calculated and demonstrates very good agreement with 
experimental data. The variation of the Raa with cen- 
trality shows a noticeable sensitivity to the presence of 
a hadronic phase. This may be used to place bounds on 
the ratio between the q in the hadronic phase compared 
to the partonic phase. More differential observables, such 
as the Raa versus the reaction plane are calculated and 
show qualitative agreement with the data at a somewhat 
lower pt- These results demonstrate the sensitivity of 
the Raa at a given angle with respect to the reaction 
plane to the initial profile of the produced matter. The 
nuclear modification factor as a function of the depth 
in the medium is shown to be anticorrelatcd with the 
length integrated density as experienced by a jet on av- 
erage. The magnitude of contributions to the Raa from 
different depths demonstrate the diminished surface bias 
in semi-central collisions. 
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